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Abstract 

High-frequency  regional  records  from  small  earthquakes  (magnitude  <  4.5),  and  comparable 
magnitude  chemical  explosions,  are  analyzed  to  find  a  reliable  seismic  discriminant  in  southern 
Russia  near  Kislovodsk.  The  digital,  three-component  seismograms  recorded  during  1992  by  the 
Caucasus  Network  operated  by  Lamont-Doherty  Earth  Observatory  since  1991  in  the  distance 
ranges  15  to  233  km  are  used.  Mean  vertical-component  PgILg  spectral  amplitude  ratios  in  the 
band  8  -  18  Hz  are  about  1.3  and  3.2  for  earthquakes  and  explosions,  respectively,  in  this  region. 
We  find  that  the  vertical-component  PgILg  spectral  ratio  in  the  frequency  band  8  - 18  Hz  serves 
quite  well  for  classifying  these  events.  A  linear  discriminant  function  analysis  indicates  that  the 
PgILg  spectral  ratio  method  provides  discrimination  power  with  a  total  misclassification 
probability  of  about  7%.  The  PgILg  spectral  ratios  of  rotated,  three-component  regional  records 
improve  the  discrimination  power  of  the  spectral  ratio  method  over  the  vertical-component  PgILg 
ratios.  Preliminary  analysis  indicates  that  distance  corrected  vertical-component  PgILg  ratios 
improve  the  discrimination  power  by  about  4%  over  uncorrected  ratios.  But  we  find  that  an  even 
better  discriminant  is  the  PgILg  spectral  ratio  of  the  three-component  regional  records  corrected 
for  the  free  surface  effect.  In  the  frequency  band  8-18  Hz,  the  free  surface  corrected  three- 
component  PgILg  spectral  ratio  provides  discrimination  power  with  a  total  misclassification 
probability  of  only  2.6%.  Free  surface  corrected  and  network  averaged  PgILg  spectral  ratios 
provide  transportability  of  the  spectral  ratio  method  to  various  regions  worldwide. 

1.  Introduction 

The  discrimination  of  small  earthquakes  from  large  chemical  explosions  (from  mines  and 
quarries)  based  on  seismic  signals  recorded  at  regional  distances  (10  -  1000  km)  is  an  important 
issue  facing  numerous  regional  seismic  networks.  The  seismic  discrimination  problem  becomes 
especially  severe  in  an  area  with  poorly  known  seismicity  such  as  southern  Russia  near 
Kislovodsk.  It  would  be  extremely  useful  to  have  a  reliable  and  robust  criterion  that  could  be  used 
to  discriminate  earthquakes  from  explosions  in  the  context  of  global  seismic  monitoring. 


2 


Figure  1  depicts  the  Caucasus  Network  operated  during  November  1991  -  summer  of  1994 
by  Lamont-Doherty  Earth  Observatory  and  Experimental  Methodological  Expedition,  Institute  of 
Physics  of  the  Earth,  Obninsk,  Russia  under  the  IRIS/Joint  Seismic  Program.  The  network 
consists  of  six  stations  around  Kislovodsk  in  southern  Russia.  The  network  was  deployed  to 
record  seismic  events  in  the  region  to  study  wave  propagation  beneath  the  Caucasus  mountains  in 
particular  and  southern  Russia  and  the  Middle  East  in  general.  The  region  covered  by  the  network 
includes  areas  that  differ  greatly  in  their  attenuation  properties,  as  previously  known  from 
WWSSN  stations  south  of  the  Caucasus  (e.g.,  Kadinsky-Cade  et  al.,  1981).  The  network  also 
has  provided  abundant  high-frequency  signals  from  small  events,  both  natural  and  man-made. 
Within  a  degree  of  the  center  of  the  network,  over  half  of  the  events  recorded  can  be  associated 
with  active  quarries,  or  roughly  150  events  per  year.  These  data  are  useful  for  characterizing 
discriminants  of  chemical  explosions  from  earthquakes,  a  critical  issue  facing  the  monitoring 
networks  planned  for  the  Comprehensive  Test  Ban  Treaty  (e.g..  Stump  et  al,  1994). 

At  each  station,  the  signals  are  recorded  on  a  three-component  set  of  intermediate-period 
sensors  with  a  natural  period,  Tq  =  5  sec  (Kinemetrics  SVl/SHl).  The  sensors  are  situated  3-4  m 
below  the  ground  surface  and  the  signals  are  digitized  on  site  by  a  digitizer  with  16  bit  resolution, 
which  after  digital  filtering  provides  a  60  samples  per  second  data  stream.  Instrument  response  of 
each  seismograph  is  nearly  flat  to  ground  velocity  in  the  frequency  band  0.2  to  24  Hz  (-3  db  level). 
Digital  signals  from  five  remote  stations  are  telemetered  to  a  central  recording  site  at  Kislovodsk 
(KIV;  Figure  1),  where  they  are  recorded  on  magnetic  tapes.  Digital  data  recorded  during  1992  are 
now  available  through  the  IRIS  Data  Management  Center. 

Most  of  the  previous  work  on  seismic  discrimination  has  focused  on  separating  large 
underground  nuclear  explosions  from  earthquakes  in  the  context  of  a  future  comprehensive  test  ban 
(e.g.,  Pomeroy  et  al,  1982;  Evernden  et  al,  1986;  Taylor  et  al,  1989).  In  most  earlier  studies, 
data  available  for  discrimination  analyses  were  limited  to  frequencies  below  10  Hz,  and  previous 
work  on  regional  signals  from  earthquakes  and  explosions  in  the  western  U.S.  suggested  that  P 
and  S  waves  from  earthquakes  often  had  stronger  high  frequencies  (that  is,  up  to  10  Hz),  than 
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signals  from  explosions  (e.g.,  Murphy  &  Bennett,  1982;  Bennett  &  Murphy,  1986;  Chael,  1988). 
There  do  not  appear  to  be  any  robust  discriminants  based  upon  regional  data  for  small  events  in  the 
passband  below  10  Hz.  However,  at  significantly  higher  frequencies,  the  PgILg  spectral  ratio 
shows  great  promise.  Thus,  Kim  et  al.  (1993)  reported  that  the  mean  PgILg  spectral  ratios  in  the 
band  1  -  25  Hz  are  about  0.5  and  1.25  for  earthquakes  and  explosions,  respectively,  in  the  eastern 
U.  S.  (A  «  10  -  610  km).  Further,  they  found  that  the  high-frequency  PgILg  spectral  amplitude 
ratio  in  the  frequency  band  5  -  25  Hz  was  an  adequate  discriminant  for  classifying  these  events.  In 
the  eastern  U.  S.,  the  PgILg  spectral  amplitude  ratio  method  provides  discrimination  power  with 
a  total  misclassification  probability  of  about  1%  (Kim  et  al.,  1993). 

The  regional  signals  from  small  seismic  events  in  the  northern  Caucasus  region  are 
characterized  by  Pg  and  Lg  onset  arrivals  with  group  velocities  of  about  5.8  km/s  and  3.2  km/s, 
respectively  at  distances  of  a  few  tens  of  km  to  about  200  km.  This  Lg  group  velocity  of  3.2 
km/s,  which  is  somewhat  slower  than  the  global  average  of  3.4  -  3.5  km/s,  indicates  lower  than 
normal  ^-wave  speeds  in  the  crust  around  the  Caucasus  Network.  Frequency  content  of  the 
regional  seismic  signals  from  the  earthquakes  and  quarry  blasts  may  be  best  examined  using 
spectrograms.  Figure  2  shows  spectrograms  and  corresponding  time  series  from  an  earthquake 
and  a  quarry  blast  (multiple-hole  instantaneous  blast)  at  a  comparable  distance  range  from  KIV  (A 
=  65  km).  The  earthquake  record  is  characterized  by  a  strong  S  wave  amplitude  at  low 
frequencies  (below  about  7  Hz),  and  P  and  S  waves  with  comparable  amplitude  at  higher 
frequencies  (10-15  Hz).  The  spectrogram  of  the  quarry  blast  record  shows  predominant  P  wave 
energy  over  a  wide  frequency  band  3  -  20  Hz  and  weak  S  waves  in  this  band  as  well  as  clear 
spectral  banding  due  to  source  multiplicity  (Figure  2).  The  time  series  of  these  events  also  show 
clear  differences  between  the  earthquake  and  explosion  records. 

Although  spectrograms  can  provide  effective  discrimination,  their  interpretation  is  usually  a 
matter  of  expert  judgement  that  cannot  be  turned  into  a  simple,  objective  algorithm  applicable  on  a 
routine  basis  to  large  numbers  of  events.  However,  recent  developments  suggest  that 
spectrograms  can  be  efficiently  used  in  seismic  discrimination  by  using  an  artificial  neural  network 
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which  can  be  trained  to  classify  earthquakes  from  explosions  (Dowla  et  al.,  1990;  Dysart  &  Pulli, 
1990  and  Wang  &  Teng,  1995).  Our  observations  of  high-frequency  (1-25  Hz)  regional  signals 
from  earthquakes  and  explosions  in  southern  Russia  indicate  that  these  seismograms  show 
distinctively  different  patterns  in  the  spectral  content  of  P  and  S  signals,  between  earthquakes  and 
explosions.  This  suggests  that  the  Pg/Lg  spectral  amplitude  ratio  can  be  used  as  a  basis  for  an 
objective  discriminant. 

In  the  following  sections,  we  report  on  our  measurements  of  the  Pg/Lg  spectral  amplitude 
ratio,  using  high-frequency  digital  seismograms  from  earthquakes  and  chemical  explosions 
recorded  by  the  Caucasus  Network  (Figure  1),  and  we  evaluate  their  discrimination  capability, 
paying  particular  attention  to  the  improvements  attainable  with  three-component  records  as 
compared  to  conventional  Pg/Lg  spectral  ratios  based  only  on  vertical-component  data. 

2.  Data 

We  present  the  analysis  of  digital  seismograms  from  25  explosions  (presumably  quarry 
blasts)  and  25  earthquakes  to  obtain  a  specific  discriminant.  Most  of  the  earthquakes  are  in  the 
magnitude  range  1  to  4  and  are  reported  by  the  Joint  Seismic  Program  Center  (JSPC)  or  the  U.S. 
Geological  Survey  (PDE).  Seismicity  during  1992,  within  a  radius  of  about  200  km  from  station 
KIV,  is  plotted  in  Figure  1.  We  examined  all  earthquakes  and  quarry  blasts  which  appear  as 
clusters  on  Figure  1 .  Events  are  selected  to  sample  a  wide  range  of  propagation  paths  around  the 
network.  All  quarry  blasts  are  from  three  known  quarry  sites.  Locations  of  mines  and  quarries  in 
the  region  are  depicted  in  Figure  3  and  listed  in  Table  1.  Accurate  locations  of  the  Tymauz  mine  in 
the  south  of  the  network  (Stump  et  al.,  1994)  as  well  as  a  Gypsum  mine  near  the  town  of  Zhako 
are  provided  by  F.  Riviere-Barbier  (Center  for  Seismic  Studies;  personal  comm.,  March  1994), 
while  locations  of  other  mines  and  quarries  are  obtained  from  a  map  provided  by  V.  Khalturin 
(personal  comm.,  March  1994).  Locations  of  the  Tyrnauz  and  Zhako  mines  obtained  from 
different  sources  as  noted  above  coincided  within  about  ±1  km  (Table  1). 

Only  explosions  for  which  we  had  information  about  the  blast  characteristics  (usually  blast 
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time  and  location;  Table  1)  or  events  whose  signals  were  clearly  identified  as  quarry  blasts  are 
included  in  the  sample  data  for  the  explosion  group.  Blast  time  is  one  of  the  most  practical 
discriminants  for  large  chemical  explosions.  About  87.5%  of  events,  77  out  of  89  events  recorded 
during  1992,  located  within  the  15  km  radius  from  the  Tymauz  mine  are  clustered  at  two  peak 
times  (10  am  &  4  pm  local  time)  while  all  3 1  events  located  within  the  10  km  radius  from  the  Ust- 
Djeguta  and  Tsementny-Zavod  quarries  are  clustered  at  2  pm  local  time  (see  Figure  1).  Distance 
ranges  of  the  data,  72  three-component  records  for  25  earthquakes  and  62  three-component  records 
for  25  explosions,  are  21  to  233  km  with  means  of  98  and  74  km  for  earthquakes  and  explosions, 
respectively.  For  vertical-component,  82  and  78  records  were  available  for  earthquakes  and 
explosions,  respectively  for  the  events. 

Pg  and  Lg  signals  are  windowed  with  a  Gaussian  weighting  function  centered  at  group 
velocities  around  5.6  km/s  and  3.2  km/s,  respectively.  A  standard  deviation,  a,  of  =  2.5  s  is 
used  for  the  Lg  Gaussian  window  at  a  reference  distance  of  Aj.^^  =100  km.  The  window  includes 
the  Lg  arrivals  with  group  velocities  between  3.48  and  2.96  km/s  in  its  ±a  width  and  the 
Gaussian  window  is  truncated  at  ±1.96  a  («  95  %),  which  include  signals  within  the  group 
velocities  3.80  to  2.77  km/s.  The  window  lengths  at  different  distances  are  scaled  by  a  =  x  A 
/  A^gf.  Thus,  the  time  window  used  for  Lg  has  the  same  group  velocity  window  regardless  of 
distance.  The  Pg  window  length  is  scaled  to  the  Lg  window  using  an  average  propagation 
velocity  ratio  VsIVp  (»  1/V3^)  x  Lg  window  length,  so  that  Pg  crustal  paths  similar  to  the  paths 
making  up  Lg  are  included. 

The  Pg  and  Lg  signals,  weighted  by  the  Gaussian  functions,  are  fast  Fourier  transformed. 
The  resulting  amplitude  spectra  are  smoothed  with  another  Gaussian  function  having  a  =  1  Hz  and 
are  re-sampled  at  every  2  Hz  interval  from  2  to  24  Hz.  Noise  analyses  indicate  that  the  signal-to- 
noise  (S/N)  ratios  are  quite  high  in  most  of  the  records,  but  in  some  cases  the  signals  fall  to  the 
background  noise  level  above  about  20  Hz.  The  S/N  ratio  becomes  less  than  2  at  about  20  Hz  for 
records  from  distant  events  (A  >  150  km).  The  log^Q(Pg/Lg)  spectral  amplitude  ratios  at  discrete 
frequency  points  are  obtained  for  each  record.  Network  averaged  log2o(Fg/Lg)  ratios  are  then 
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obtained  for  each  event  by  averaging  the  discrete  frequency  values  from  all  stations.  Figure  4 
shows  the  network  averaged  logjQ(Pg/Lg)  ratios  at  discrete  frequency  points  for  25  earthquakes 
and  25  explosions.  Earthquake  and  explosion  populations  are  well  separated  in  the  frequency  band 
6  to  18  Hz  (Figure  4).  At  frequencies  lower  than  6  Hz,  there  is  some  overlap.  Thus,  the 
frequency  band  6  -  18  Hz  is  the  most  reliable  and  a  sampled  version  will  be  used  in  the  following 
discrimination  analysis. 

3.  The  Linear  Discriminant  Function 

To  test  the  discriminant  power  of  the  high-frequency  PgILg  spectral  amplitude  ratio,  we 
performed  multivariate  discriminant  analysis  on  vertical-component  logjQ(Pg/Lg)  spectral  ratios 
for  the  data  set  of  earthquakes  and  explosions.  The  discriminant  power  is  evaluated  by  developing 
an  optimal  classification  rule  that  separates  earthquakes  from  explosions  based  on  the  sample  data 
we  have  on  each  group.  Each  sample  data  set  of  25  events  (i.e.  25  explosions,  25  earthquakes)  is 
described  by  a  matrix  of  6  rows  [(logjQ(Pg/Lg)  at  2  Hz  frequency  intervals  from  8  to  18  Hz]  and 
25  columns.  Preliminary  analysis  suggested  somewhat  higher  PgILg  ratios  from  Tymauz  mine 
blasts  than  for  other  quarry  blasts,  but  we  merged  both  data  sets  into  one  explosion  group  because 
the  differences  between  them  were  much  smaller  than  their  differences  from  earthquakes. 

We  introduce /£;^(r)  and/^^Cr)  as  the  probability  densities  of  the  two  types  of  events,  with  r 
as  a  column  vector  representing  the  log  spectral  ratios  sampled  at  6  values.  We  take  and 
as  the  a  priori  probability  of  the  two  types  of  events,  so  +  ^Eq  -  1  •  We  follow  standard 
practice  by  assigning  an  event  to  the  earthquake  class  if 

fEqi^)  I  fExi^)  >  ^Ex  I  ^Eq  (1) 

and  to  the  explosion  class  otherwise  (see  e.g.,  Seber,  1984  for  multivariate  discrimination 
methods).  This  rule  is  optimum,  in  that  it  minimizes  the  total  probability  of  misclassification.  Our 
knowledge  of^^^fr)  and/£^(r)  comes  from  the  sample  data  set,  and  we  introduce  and  evaluate  a 
linear  discriminant  function  D(r)  under  the  assumptions  that:  a)  the  sample  distributions  are 
normal;  b)  the  dispersion  (variance-covariance)  matrices  of  the  two  groups  are  the  same;  and  c)  the 


7 


sample  data  observations  are  correctly  classified.  Then  the  linear  function  that  best  separates  the 
two  types  of  seismic  source  is 

Z)(r)=A''[r-(H^^+MgP/2]  (2) 

where  is  the  mean  of  the  earthquake  sample  data  values  {i  =  1,...,25),  and  similarly  for 
;  A  =  S“'  -  Pg.^ ),  and  S  is  the  average  of  the  dispersion  matrices  and 

1  25  .  T  .  T 

with  (for  example)  ^  Z  (JEq^  ~  denotes  a  transpose.  Equation 

(2)  defines  a  hyperplane  for  separating  the  two  groups. 

In  terms  of  these  easily  calculated  quantities,  the  discrimination  analysis  is  very  simple.  The 
rule  (1)  becomes:  assign  an  event  r  to  the  earthquake  population  if  D(r)  >  In  /  tce^).  If 
=  Kex  =  j  earthquakes  and  explosions  have  equal  a  priori  probabilities)  the  rule  is  even 
simpler:  the  event  is  labelled  an  earthquake  or  an  explosion  according  to  D{r)  >  0  or  D(r)  <  0.  In 
this  case,  we  have 

,  r-A/2  0  T 

P  (misclassification)  =  exp(-x2/2)  dx  where  A  =  A  (p^.^  -  Pg,^ )  (3) 

for  the  probability  of  misclassification.  is  the  Mahalanobis  D-squared  measure  of  distance 
between  the  two  means. 

This  type  of  multivariate  analysis  was  developed  and  used  by  Fisher  (1936),  and  has  been 
used  in  many  seismological  applications  —  for  example  by  Jebe  and  Willis  (1964),  Sandvin  and 
Tjpstheim  (1978)  and  Kim  et  al.  (1993).  The  equality  of  the  dispersion  matrices  for  the  two 
groups  is  tested  using  a  test  (Seber,  1984).  If  the  dispersion  matrices  of  the  populations  are  not 
the  same,  then  a  quadratic  discriminant  function  is  appropriate  (Davis,  1986;  Seber,  1984). 

4.  Applications  of  Discriminant  Analysis 

4.1  High-Frequency  Network  Averaged  Vertical-Component  PgILg  Ratio 

The  sample  data  sets  consisting  of  25  earthquakes  and  25  explosions  were  analyzed  using  the 
linear  discriminant  function  given  in  (2).  For  each  event,  network  averaged  vertical-component 
logjQ(Pg/Lg)  ratios  at  frequencies  of  8,  10,  12,  14,  16  and  18  Hz  correspond  to  the  variables  r^ 
^2’  ^3’  ^4’  ^5  ^6'  linear  discriminant  function  obtained  is 
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D(r)  =  6.68  +  10.04  -  S.22r^  -  2.03r^  -  IS.SSr^  +  6.l4r^  - 1  l.lSr^  (4) 

and  the  Mahalanobis  D-squared  measure  is  =  8.67.  Assuming  equal  prior  probabilities  for  the 
two  groups,  we  assign  event  r  to  the  earthquake  class  if  D(r)  >  0.  Applying  this  rule  to  the 
earthquake  and  explosion  data,  we  find  that  three  earthquakes  and  an  explosion  are  classified 
incorrectly  and  the  misclassification  probability  from  eq.  (3)  is  0.071.  Values  of  D(r)  may  be 
called  the  discriminant  score  and  are  plotted  in  Figure  5  with  respect  to  the  mean  logjQ(Pg/Lg) 
spectral  amplitude  ratio  of  each  event.  Vertical  lines  in  the  figure  denoted  as  Eq  and  Ex  are  the 
projection  of  the  multivariate  mean  of  the  earthquake  and  explosion  populations,  respectively.  The 
vertical  line,  Dq,  is  the  line  D(r)  =  0,  which  serves  to  classify  the  events  when  the  a  priori 
probability  of  the  two  populations  is  the  same.  The  distance  between  Eq  and  Ex  is  the 
Mahalanobis  D-squared  measure  of  distance  between  two  populations,  since,  from  (2),  - 

D(\Iex)  =  A^.  It  is  shown  in  Figure  5  that  all  the  earthquake  records  from  various  paths  in  the 
southern  Russia  near  Kislovodsk  have  a  mean  vertical-component  Pg/Lg  spectral  ratio  of  about 
1.3,  while  the  explosion  records  show  a  mean  of  about  3.2  in  the  frequency  bands  8-18  Hz.  The 
locations  of  three  misclassified  earthquakes,  -  events  #8,  #1 1  and  #16  (see  Table  2)  -  are  indicated 
in  Figure  5. 

4.2  Network  Averaged  Three-Component  Pg/Lg  Ratio 

The  Caucasus  Network  and  many  other  modem  seismograph  stations  are  now  equipped  with 
three-component  broadband  sensors.  Although  high-frequency  regional  signals  on  the  vertical- 
component  show  reasonable  discrimination  power  for  classifying  earthquakes  from  explosions, 
there  is  interest  in  the  evaluation  of  three-component  records  to  see  if  they  can  provide  improved 
discrimination.  We  have  therefore  examined  Pg/Lg  spectral  ratio  methods  to  discriminate  the 
earthquakes  from  explosions  using  regional  three-component  records.  Our  first  approach  was 
simply  to  rotate  the  observed  NS-  and  EW-component  seismogram  pairs  to  obtain  radial  (R)  and 
tangential  (T)  components.  Then  the  Pg/Lg  ratios  of  three-component  records  are  formed  for 
each  station  by  defining 
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<jLgl+Lgl+Lgl 

where  subscripts  indicate  the  component.  For  example,  Pg  ^  indicates  spectral  amplitude  of  P 
waves  (mostly  Pg)  on  vertical-component  (Z),  Lg  ^  indicates  spectral  amplitude  of  Lg  waves  on 
radial-component.  A  single  three-component  Pg/Lg  ratio  is  obtained  for  each  station  from  eq.  (5) 
and  the  network  mean  logjQ(Fg/Lg)  is  obtained  by  averaging  over  stations  for  each  event.  Figure 
6  shows  the  network  averaged  logjQCPg/Lg)  ratios  at  discrete  frequency  points  for  25 
earthquakes  and  25  explosions.  Earthquake  and  explosion  populations  are  better  separated  in  the 
frequency  bands  2-4  and  12-20  Hz  than  the  vertical-component  ratios  shown  in  Figure  4.  There  is 
also  a  clear  difference  between  the  spectral  ratios  at  low  frequencies  (2-10  Hz)  and  high 
frequencies  (12-20  Hz)  within  the  same  group,  in  particular,  for  earthquake  data  (Figure  6).  The 
network  averaged  log^Q(Pg/Lg)  ratios  of  the  three-component  records  from  the  earthquakes  along 
various  paths  in  southern  Russia  have  a  mean  three-component  Pg/Lg  spectral  ratio  of  about 
0.73,  while  the  explosion  records  show  a  mean  of  about  2.2.  These  mean  Pg/Lg  ratios  are 
slightly  lower,  by  about  a  factor  of  1.6,  than  the  vertical-component  data  shown  earlier  (see  Figure 
4),  due  mainly  to  the  inclusion  of  Lg  waves  on  radial  and  tangential  components  (Lg  is  typically 
stronger  on  these  components,  than  on  the  vertical-component). 

We  repeated  the  multivariate  analysis  using  three-component  Pg/Lg  ratios.  The  discriminant 
scores  of  each  event  are  plotted  with  respect  to  the  network  mean  logjQ(Pg/Lg)  spectral  amplitude 
ratio  in  Figure  7  and  are  listed  in  Table  2  and  Table  3.  Discriminant  function  analyses  for  the  three- 
component  data  suggest  that  the  use  of  three-component  records  reduces  the  scatter  of  the  mean 
logjQ(Pg/Lg).  The  Mahalanobis  D-squared  measure  is  =  14.07  indicating  a  greater  separation 
between  earthquake  and  explosion  populations  than  was  the  case  with  vertical-component  data. 
Three  distant  earthquakes,  events  #8,  #1 1  and  #16  (see  Figure  3),  that  were  misclassified  in  the 
vertical-component  analysis  are  now  correctly  classified.  In  general,  distant  earthquakes  show 
higher  mean  log2Q(Pg/Lg)  than  near  events  (Figure  3  and  Figure  7).  Most  of  the  Tyrnauz  mine 
explosions  have  higher  mean  spectral  ratios  than  signals  from  the  other  two  mining  areas  (Figure 
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7)- 

The  discriminant  function  analysis  indicates  that  all  events  in  the  sample  data  are  classified 
correctly  and  the  total  misclassification  probability  is  only  3.04%,  suggesting  that  the  three- 
component  PgILg  ratios  improve  the  discrimination  power  of  high-frequency  regional  data  by 
about  4%  over  the  use  of  single  vertical-component  Pg/Lg  ratios  (92.9%  versus  97.0%). 

5.  Removal  of  Free  Surface  Effect  on  Regional  Records 

Although  high-frequency  Pg/Lg  spectral  ratios  of  regional  signals  on  vertical-component 
records  show  reasonable  discrimination  power  for  classifying  earthquakes  from  explosions  -  and 
rotated  three-component  records  are  apparently  a  significant  improvement  -  the  success  of  such 
discriminants  in  trial  regions  does  not  always  warrant  their  applicability  to  various  other  regions 
world- wide.  The  key  to  thie  successfiil  discrimination  of  various  types  of  seismic  sources  is  the 
clear  observation  of  the  signals  radiated  from  the  seismic  sources.  One  of  the  most  important 
problems  is  correction  of  the  effects  of  the  source-receiver  paths  on  observed  regional  signals.  The 
source-receiver  path  effects  include,  for  example,  attenuation  along  the  paths,  different  geometrical 
spreading  of  P  and  S  waves,  local  receiver  site  responses  and  interaction  of  the  incoming  wave 
field  with  the  free  surface. 

The  effect  of  the  free  surface  in  the  neighborhood  of  the  receiver  on  the  recorded  regional 
signals  is  significant  and  its  main  effects  are  amplification  of  incoming  waves.  However,  unlike 
other  path  effects,  it  can  be  removed  by  a  relatively  simple  real  transformation  (House  & 
Boatwright,  1980;  and  Kennett,  1991).  Recently,  Kennett  (1991)  showed  that  a  set  of 
approximate  free  surface  correction  operators  can  be  formed  to  remove  the  free  surface  effects  over 
the  slowness  bands  for  the  main  regional  phases. 

For  a  P-wave,  with  pulse  shape  P(t)  and  particle  motion  having  (jc,  y,  z)  components 
P(0  (sin  /,  0,  -cos  i),  incident  upon  a  free  surface  perpendicular  to  the  z  (depth)  direction,  the 
motion  of  the  free  surface  itself  is  (Aki  and  Richards,  1980;  page  190), 
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where  a  =  surface  P-wave  velocity,  =  surface  i'-wave  velocity,  i  =  angle  of  incidence  of  the 

P-wave, ;  =  angle  of  incidence  of  the  5-wave,  p  =  slowness  and  sin  Ua  =  sin  jlfi  =  p. 

Note  that  Z(t)  is  conventionally  measured  positive  upward,  opposite  to  the  depth  direction;  and 

P(0  is  positive  away  from  the  source.  The  ground  motions  in  a  homogeneous  half-space  for  a 

P-wave  incident  upon  a  free  surface  are  depicted  in  Figure  8. 

The  motion  of  the  incident  P-  and  5V-wave  can  then  be  recovered  as  a  linear  combination 

of  the  radial,  Rit),  and  vertical,  Z(0,  components  available  from  three-component  records. 

Thus,  equation  (6)  can  be  used  to  show  that 
cos  2  7  B 

Pit)  =  - - l.ZfO  +  ^sin;  .P(0.  (7a) 

2  cos  I  « 

The  corresponding  result  for  an  incident  5V-wave  is 

SVit)  =  ^^-Rit)-smj  ^Z(t),  (7b) 

2  cos ; 

which  has  long  been  used  by  those  studying  earthquake  source  functions  to  infer  the  5-wave 
pulse  shape  from  three-component  data  (House  and  Boatwright,  1980). 

The  incident  5/f-wave  is  simply  one  half  of  the  recorded  tangential-component: 

5/7(0  =  ^  T(0.  (7c) 

If  the  incident  P-wave  is  contaminated  by  SV  motion  due  to  local  horizontal  layering,  then 
eq  (7a)  gives  precisely  the  right  combination  of  R  and  Z  to  eliminate  the  contribution  from  5V. 
Similarly,  eq.  (7b)  removes  contaminating  P  (due  to  horizontal  layering)  from  the  SV  trace. 

An  example  of  three-component  records  of  an  explosion  from  the  Tyrnauz  mine  recorded  at 
station  MIC  (A  =  108  km)  is  shown  in  Figure  9.  The  rotated  seismograms,  Z,  R  and  T 
components,  show  that  P  waves  on  the  Z-component  as  well  as  Lg  waves  on  the  R-  and  T- 
components  have  comparable  amplitudes  (Figure  9).  However,  when  free  surface  effects  are 
removed  from  the  rotated  seismograrhs  by  successively  applying  free  surface  correction  operators 
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appropriate  for  major  regional  phases  (Pn,  Pg  and  Lg),  the  recovered  incident  wavevector 
estimates  allow  us  directly  to  compare  the  P-,  SV-  and  SH-wave  amplitudes  on  different  traces. 

The  composite  incident  wavevector  traces  in  Figure  9  are  obtained  by  applying 
transformations  (equations  7abc)  sequentially  to  successive  group  velocity  windows  along  the 
trace.  Successive  windows  used  are  indicated  in  Figure  9  with  corresponding  slownesses. 
Assuming  a  =  4.5  km/s  and  =  2.6  km/s  appropriate  for  the  region,  signals  with  group 
velocities  greater  than  5.2  km/s  are  transformed  with  an  approximate  Pn  slowness,  p  =  0.08  s/km 
(/  «  25°),  signals  with  group  velocities  5.2  -  4  km/s  are  transformed  with  an  approximate  Pg 
slowness,  p  =0.14  s/km  (/  «  40°)  and  signals  with  group  velocities  3.3  -  2.6  km/s  are 
transformed  with  an  approximate  Lg  phase  slowness,  p  =  0.34  s/km  (/  «  60°).  The  SV 
transformation  is  sensitive  to  the  slowness  used,  since  the  free  surface  correction  factor  for  R- 
component  in  (7b)  varies  rapidly  for  slownesses  l/(^/2j8)  <  p  <  l/p.  When  the  precise 
slownesses  of  the  incoming  waves  are  not  available  for  correct  transformation,  our  experience 
suggests  that  a  reasonable  choice  of  the  slowness  is  when  the  P-wave  on  the  SV-trace  after  the 
transformation  has  comparable  amplitude  to  the  P-wave  on  the  5/7-trace,  since  P  waves  on  both 

traces  are  derived  from  scattering,  for  example  from  heterogenieties  other  than  horizontal  layering. 
Similarly  for  Lg  signals  from  the  explosions,  a  reasonable  choice  of  the  slowness  is  when  the  Lg 

waves  on  the  P-trace  are  comparable  to  the  Lg  waves  on  the  5/7-trace.  For  p  >  1/a,  the 

coefficient,  i§  complex,  since  P-waves  become  evanescent  and  cos  j  =  (1  - 

2  cos  i 

is  imaginary.  In  this  slowness  range  there  is  a  phase  shift  introduced  between  the  vertical  and 
radial  components  of  the  5^  contributions  to  the  surface  displacement  which  gives  rise  to  elliptical 
polarization  patterns  (Kennett,  1991).  Therefore,  for  the  P-trace,  the  Z-component  is  Hilbert 
transformed  before  the  transformation  (7a)  for  slownesses  greater  than  1/a. 

The  composite  incident  wavevector  traces  show  that  indeed  the  P  waves  are  mainly  on  the 
P-wave  vector  component,  while  5  waves  (Lg)  are  dominantly  on  the  5  V- wave  vector 
component  (Figure  9),  hence  the  amplitudes  of  major  crustal  phases  reveal  clearer  radiation 
characteristics  from  the  source.  In  particular,  the  relatively  small  5/7  in  the  regional  5  waves 
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clearly  indicates  the  signal  being  generated  by  an  explosion  source.  A  very  useful  result  of  the  free 
surface  correction  is  that  we  can  make  direct  comparison  between  the  P-,  SV-  and  SH-wave 
amplitude  in  a  particular  group  velocity  window  and  so  get  closer  to  the  radiation  characteristics 
from  the  source.  This  is  of  potential  significance  for  discriminating  different  types  of  sources  and 
we  examine  the  usefulness  of  free  surface  corrected  data  in  the  following  section. 

5.1  Network  Averaged  Free  Surface  Corrected  Three-Component  Pg/Lg  Ratio 

We  examined  the  PgILg  spectral  ratio  method  to  discriminate  earthquakes  from  explosions 

using  the  free  surface  corrected  regional  three-component  records.  The  PgILg  ratios  from  the 

free  surface  corrected  P-  and  SV-  and  5H-seismogram  of  three-component  regional  records  are 

obtained  for  each  station  by  defining 

Ps  _  PSp 
Lg  jjl  \  1 

where  Pg  p  =  spectral  amplitude  of  the  P-wave  on  the  P-seismogram,  Lg  =  spectral  amplitude 
of  the  Lg-wave  on  the  5  V-seismogram  and  Lg  =  spectral  amplitude  of  the  Lg-wave  on  the 
5H-seismogram.  A  single  PgILg  ratio  is  obtained  for  each  station,  and  the  network  mean 
PgILg  ratio  is  obtained  by  averaging  over  stations  for  each  event.  The  network  averaged 
logjQ(Pg/Lg)  ratios  based  on  eq.  (8)  from  P-,  SV-  and  5/f-seismograms  are  plotted  in  Figure 
10.  The  separation  of  the  spectral  ratios  of  earthquake  and  explosion  groups  appears  to  be  similar 
to  the  three-component  Pg/Lg  ratios  (Figure  6). 

The  sample  data  sets  consisting  of  25  earthquakes  and  25  explosions  were  analyzed  using  the 
linear  discriminant  function  as  in  the  previous  sections.  For  each  event,  network  averaged  free 
surface  corrected  three-component  log,Q(Pg/Lg)  ratios  at  frequencies  of  8,  10,  12,  14,  16,  18 
and  20  Hz  correspond  to  the  variables  r^,  r2,  r^  r^,  r^  and  r^.  The  linear  discriminant 
function  that  best  separates  the  two  types  of  source  is 

£>(r)  =  -5.09  +  7.17rj  -  17.43r2  +  17.52r3  -  32.52r4  +  11.59^5  -  \2.l5r^  (9) 

and  the  Mahalanobis  D-squared  measure  of  distance  between  the  two  population  means  is  = 
15.15,  higher  than  either  of  the  two  values  obtained  above.  Assuming  equal  prior  probabilities  for 
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the  two  groups,  and  applying  the  previous  classification  rule  to  the  earthquake  and  explosion  data, 
we  find  that  two  explosions  are  classified  incorrectly  and  the  total  misclassification  probability  is 
very  low,  namely  2.6%.  The  two  misclassified  explosions  are  events  #21  and  #25  (Table  3). 

The  discriminant  score  for  each  event  is  plotted  in  Figure  1 1  with  respect  to  the  network 
mean  free  surface  corrected  log^Q{Pg/Lg)  spectral  amplitude  ratio  of  each  event.  It  is  shown  in 
Figure  1 1  that  the  earthquake  records  from  various  paths  in  southern  Russia  have  a  mean  free 
surface  corrected  three-component  Pg/Lg  spectral  ratio  of  about  0.32  in  the  frequencies  8-18  Hz, 
while  the  explosion  records  show  a  mean  of  about  1.0.  These  mean  Pg/Lg  ratios  are  lower  than 
the  vertical-component  or  three-component  data  discussed  earlier  (see  e.g..  Figure  5  &  7)  by  about 
a  factor  2  to  3  and  there  is  a  clearer  separation  of  Pg/Lg  spectral  ratios  of  the  earthquake  and 
explosion  populations.  Notice  that  two  distant  earthquakes,  #2  &  #16,  have  higher  Pg/Lg  ratios 
than  other  earthquakes  (Figure  11).  The  free  surface  corrected  Pg/Lg  spectral  ratios  at 
frequencies  5  -  20  Hz  indicate  that  for  the  earthquakes,  the  radiated  seismic  amplitudes  are 
partitioned  into  P  and  SV-SH  waves  with  a  ratio  1:3,  while  for  the  explosions,  P  and  SV  waves 
have  comparable  amplitudes.  Use  of  the  free  surface  corrected  Pg/Lg  spectral  ratio  improves  the 
discrimination  power  of  high-frequency  regional  data  by  about  4.5%  over  the  use  of  single 
vertical-component  ratios  (92.9%  versus  97.4%). 

6.  Correction  of  Distance  Effects  on  the  Pg/Lg  Spectral  Amplitude  Ratio 

The  sample  data  -  25  earthquakes  and  25  explosions  -  used  for  the  linear  discriminant 
function  analyses  in  the  previous  sections  come  from  unevenly  distributed  sources.  The  explosions 
were  exclusively  from  three  mining  areas,  while  the  earthquakes  were  widely  distributed  in  the 
region  (Figure  3).  For  simplicity,  we  could  have  chosen  only  earthquakes  that  were  located  at 
distance  ranges  comparable  to  the  explosions,  however  we  opted  to  select  earthquakes  at  various 
locations  to  sample  a  representative  seismicity  of  the  region,  for  developing  seismic  discriminants 
for  practical  monitoring  purposes  in  the  region.  Because  the  distribution  of  the  earthquakes  and 
explosions  used  were  uneven,  we  examined  the  effects  of  distance  on  the  spectral  ratio. 
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The  vertical-component  PgILg  ratios  of  82  earthquake  records  and  78  explosion  records  are 
plotted  against  distance  in  Figure  12.  The  mean  of  the  PgILg  ratios  in  the  frequency  band  8-18 
Hz  are  plotted  and  an  average  trend,  calculated  from  each  of  the  distance  bins  in  every  20  km 
interval,  is  superposed  in  the  figure.  The  earthquake  data  show  a  weak  distance  dependence  in  the 
ratio,  caused  probably  by  differences  in  attenuation  of  Pg  (or  P)  and  Lg  (or  S)  phases  over 
distance.  The  PgILg  ratios  at  distances  less  than  about  100  km  are  about  -0.1  logarithmic  units 
(l.u.)  and  they  increase  to  0.2  l.u.  at  distance  ranges  of  100  to  180  km.  The  ratios  reach  to  0.4  l.u. 
at  distances  180-240  km  (Figure  12).  Hence,  the  PgILg  ratios  of  earthquake  data  increase  by  a 
factor  of  about  three  within  the  distance  range  of  80  to  200  km.  We  notice  a  spectral  ratio 
minimum  in  the  earthquake  data  at  about  170  km.  It  is  probably  associated  with  the  crossover 
distance  in  P  wave  travel  time  curve  where  Pn  phase  becomes  the  first  arrival  phase.  It 
represents  a  change  in  the  regional  wave  field  in  which  P  waves  start  leaking  to  the  upper  mantle, 
and  thereby  reducing  P  eiiergy  compared  to  S  or  Lg.  This  spectral  ratio  minimum  occurs  at 
approximately  180  km  for  the  northeastern  U.  S.  data  set  studied  by  Kim  et  al.  (1993).  It  is  noted 
that  the  “distance”  here  does  not  represent  a  simple  length  scale,  rather  it  is  the  “path  length”  within 
a  complex  3-D  earth  structure,  since  the  paths  from  most  of  the  distant  earthquakes  cross  a  major 
tectonic  boundary  -  Greater  Caucasus  (see  Figure  1  &  3).  Nevertheless,  we  examine  a  first-order 

approximation  to  account  for  the  effects  of  “distance”  on  the  spectral  ratio. 

For  Pg  and  Lg  phases,  the  spectral  amplitude  decay  with  distance  can  be  written  in  the 

form, 

A(A,/)  =  Ao(/)  exp  (-Y(/)  A)  A-^  (10) 

which  consists  of  an  anelastic  attenuation  factor  and  the  geometrical  spreading  factor, where  A  = 
epicentral  distance  in  km,  A(A,/)  =  observed  spectrum  at  distance  A,  A^if)  =  source  spectrum, 
y{f)  =  anelastic  attenuation  coefficient  and  k  =  exponent  for  geometrical  spreading.  7(/)  is  related 
to  quality  factor  Qif)  by  the  formula:  yif)  =  itf  l{Q(f)U),  where  U  is  the  group  velocity. 
Q(f)  is  a  function  of  frequency,/,  and  is  assumed  to  have  the  form:  Q(f)  =!2o/^»  where  Qq 
is  the  quality  factor  at  a  reference  frequency  (f=  1  Hz)  and  is  the  exponent  of  frequency.  Using 
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the  differences,  y(f)  =  yp(f)  -  YsC/)  and  k  =  Kp  -  the  distance  corrected  spectral  ratio, 
^ogiQ{Pgo(f)/Lgo(f)},  can  be  written  as 

^ogio{Pgo(f)/LgQ(f)}  =  logjo  [{Pg(f,  ^)/Lg(f,  A)}  exp  (y(f)  A)  A’'] .  (11) 

A  recent  study  by  Abers  &  Sarker  (1996)  for  the  northern  Caucasus  region  suggests  average 
frequency  independent  Q  values  of  10701200  and  21301600  for  Pg  (or  P)  waves  in  the 
predominantly  mountain  and  shield  path,  respectively.  They  obtained  Q  values  of  775175  and 
20601315  for  Lg  (or  S)  waves  in  the  mountain  and  shield  path,  respectively  in  the  frequencies  up 
to  10  Hz.  The  boundary  between  the  two  tectonic  features  (i.e.,  mountain  and  shield)  lies 
approximately  at  the  northern  edge  of  the  Greater  Caucasus  (Figure  1)  and  it  closely  follows  the 
northern  edge  of  the  2000  m  contour  line  shown  in  Figure  3.  These  Q  values  for  Lg  waves  with 
U  =  3.4  km/s  give  y(1  Hz)  =  0.0012  km“l  and  0.00045  km“l  for  mountain  and  shield  path, 
respectively.  These  attenuation  coefficients  seem  to  be  somewhat  small  when  compared  with  other 
more  stable  regions.  For  instance,  a  recent  study  by  Shi  et  al.  (1996)  indicates  an  average  Qq  « 
625  and  r\  ~  0.25  for  Lg  waves  in  the  northeastern  U.  S.,  which  give  y(1  Hz  )  =  0.00148  km“l  at 
a  reference  frequency  of  1  Hz  with  U  =  3.4  km/s.  For  Pg  waves,  above  Q  values  with  U  =  6 
km/s  give  Y(1  Hz)  =  0.00049  km“i  and  0.00025  km“i  for  mountain  and  shield  path,  respectively. 

We  assume  that  the  geometrical  spreading  for  both  P  and  S  waves  at  short  distances  (A  < 
100  km)  is  similar  and  that  it  may  be  approximated  by  A“k  The  geometrical  spreading  for  spectral 
amplitude  of  Lg  waves  ranges  from  A“®-^  to  A~1  (e.g.,  Campillo  et  al.,  1983),  while  Pg  waves 
range  between  A“l  and  A~l-5  (e.g.,  Langston,  1983)  at  regional  distances.  These  values  indicate  a 
differential  geometrical  spreading  factor  with  an  exponent  of  ic  =  0  for  A  <  100  km  and  ic  «  0.5 
for  A  >  100  km. 

Combining  the  above  differential  geometrical  spreading  and  attenuation  coefficient  between 
Pg  and  Lg  waves,  the  distance  correction  for  the  spectral  ratio  in  northern  Caucasus  region  can  be 
written  as  exp(-0.0002  knr^fA)  and  exp(-0.0007  km“^/A)  A*^  for  shield  and  mountain 
path,  respectively.  These  distance  correction  factors  are  plotted  in  Figure  13  against  distance.  We 
indicate  as  “mixed  path”  the  propagation  paths  consisting  of  shield  and  mountain  paths  in  about 
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equal  parts. 

These  distance  corrections  are  applied  to  the  vertical-component  spectral  ratios  for  each  path, 
and  resulting  ratios  are  plotted  in  Figure  12  against  epicentral  distance  for  comparison.  As 
expected  from  Figure  13,  the  corrected  ratios  in  the  high-frequency  band  (8-18  Hz)  are  about  0.1 
l.u.  lower  for  earthquake  records  at  a  distance  of  100  km,  while  the  corrected  ratios  are  about  0.4 
l.u.lower  than  the  raw  uncorrected  ratios  at  a  distance  of  240  km  (Figure  12).  The  trend  after  the 
correction  shows  a  much  weak  distance  dependence  of  the  ratios  and  suggests  that  the  distance 
correction  using  eq.  (11)  may  be  a  reasonable  first-order  approximation  to  account  for  distance 
effects  on  the  spectral  ratio. 

6.1  Distance  Corrected  Network  Averaged  Vertical-Component  PgILg  Ratio 

We  examined  the  PgILg  spectral  ratio  method  to  discriminate  earthquakes  from  explosions 
using  distance  corrected  regional  vertical-component  records.  The  PgILg  ratios  from  vertical- 
component  records  corrected  for  distance  by  using  eq.  (11)  are  used  as  the  sample  data  set.  A 
single  PgILg  ratio  is  obtained  for  each  station,  and  the  network  mean  PgILg  ratio  is  obtained  by 
averaging  over  stations  for  each  event.  The  network  averaged  logy^lPglLg)  ratios  corrected  for 
the  distance  are  plotted  in  Figure  14.  The  mean  and  standard  deviation  at  each  frequency  is 
superposed  in  the  figure  to  indicate  the  quality  of  the  distance  corrected  ratios.  The  distance 
corrected  spectral  ratios  show  smaller  PgILg  values  and  smaller  scatter  than  the  uncorrected  ratios 
at  all  frequencies  up  to  about  20  Hz.  The  separation  of  the  spectral  ratios  of  earthquake  and 
explosion  groups  appears  to  be  similar  to  the  uncorrected  ratios  at  low  frequencies  (<  8  Hz),  but 
the  separation  between  the  two  groups  is  increasing  with  frequency  for  frequencies  greater  than  10 
Hz.  This  is  due  to  the  fact  that  several  earthquakes  were  more  distant  than  explosions. 

For  frequencies  higher  than  20  Hz,  corrected  ratios  show  greater  scatter  than  the  raw  ratios, 
suggesting  that  the  distance  correction  applied  was  not  appropriate  at  these  high  frequencies 
(Figure  14).  This  is  probably  due  to  fact  that  the  Q  factors  reported  by  Abers  &  Sarker  (1996)  are 
obtained  for  frequencies  up  to  10  Hz.  The  frequency  independent  Q  factors  we  have  used  may 
not  be  applicable  at  frequencies  higher  than  20  Hz  as  shown  by  Shi  et  al.  (1996)  who  showed  the 
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necessity  for  frequency  dependent  Q  in  regional  waves  studied  up  to  35  Hz.  However,  we  expect 
the  frequency  independent  Q  we  have  used  in  the  present  study  provides  fairly  reliable  spectral 
ratios  for  frequencies  up  to  18  Hz  as  shown  in  Figure  14  and  a  sampled  version  will  be  used  in  the 
following  discrimination  analysis. 

The  sample  data  sets  consisting  of  25  earthquakes  and  25  explosions  (vertical-component 
only)  were  analyzed  using  the  linear  discriminant  function  as  in  the  previous  sections.  For  each 
event,  distance  corrected  network  averaged  logjQ(Pg/Lg)  ratios  at  discrete  frequencies,  8-18  Hz, 
are  used  as  before.  The  linear  discriminant  function  that  best  separates  the  two  types  of  source  is 

D(r)  =  6.06  +  16.67rj  -  15.57r2  -  9A6r^  -  18.33r4  +  A.'i6r^  -  6.93r^  (12) 

and  the  Mahalanobis  D-squared  measure  of  distance  between  the  two  population  means  is  = 
14.14,  much  higher  than  the  value  obtained  for  distance  uncorrected  case  (cf.  A^=  8.67). 
Assuming  equal  prior  probabilities  for  the  two  groups,  and  applying  the  previous  classification  rule 
to  the  earthquake  and  explosion  data,  we  find  that  an  earthquake  is  classified  incorrectly  and  the 
total  misclassification  probability  is  3%.  This  is  very  close  to  the  misclassification  probability 
obtained  for  the  three-component  analysis.  Records  from  a  single  misclassified  earthquake  (event 
#16)  show  much  stronger  Pg  than  Lg  signals  on  vertical-components.  As  such  the  event  could 
only  be  correctly  classified  by  using  three-component  data. 

By  examining  a  variety  of  different  discrete  frequency  bands  using  linear  discriminant 
function  analysis,  we  found  that  the  frequency  band  6-14  Hz  provided  the  best  discrimination 
power  for  this  data.  Apparently  the  distance  correction  has  increased  the  discrimination  power  for 
the  frequency  band  8-14  Hz,  while  reducing  the  discrimination  power  for  higher  frequency  band 
16-18  Hz. 

The  discriminant  score  for  each  event  is  plotted  in  Figure  15  with  respect  to  the  network 
mean  distance  corrected  logjQ(Pg/Lg)  spectral  ratio  of  each  event.  The  mean  Pg/Lg  ratios  are 
lower  than  the  distance  uncorrected  vertical-component  (see  e.g..  Figure  5)  and  there  is  a  clearer 
separation  of  Pg/Lg  spectral  ratios  of  the  earthquake  and  explosion  populations  than  uncorrected 
case  (cf.  Figures  5  &  15).  This  preliminary  analysis  indicates  that  use  of  the  distance  corrected 
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PgILg  spectral  ratio  improves  the  discrimination  power  of  high-frequency  regional  data  by  about 
4.1%  over  the  use  of  single  vertical-component  ratios  (92.9%  versus  97.0%).  However,  more 
accurate  information  regarding  the  geometrical  spreading  and  anelastic  attenuation  of  the  regional 
phases  such  as  Pg  and  Lg  are  required  for  applying  the  distance  correction  to  PgILg  ratios  with 
confidence. 

1.  Discussion  and  Conclusions 

We  find  that  Pg  waves  and  Lg  waves  from  earthquakes  have  comparable  high  frequency 
content  over  the  broad,  high-frequency  band  width  of  5  -  20  Hz  at  regional  distances  in  southern 
Russia  near  Kislovodsk:  the  mean  vertical-component  PgILg  ratio  for  82  earthquake  records  is 
about  1.3.  For  explosions  (multiple-hole  instantaneous  and  ripple-fired  explosions),  the  Pg 
waves  are  much  stronger  than  Lg:  the  mean  vertical-component  PgILg  ratio  for  78  explosion 
records  in  the  same  frequency  band  is  about  3.2.  These  mean  PgILg  ratios,  1.3  for  earthquakes 
and  3.2  for  explosions,  are  much  higher  than  the  ratios  observed  in  tectonically  stable  eastern  U. 
S.,  where  the  ratios  are  0.5  and  1.25  for  the  earthquakes  and  explosions,  respectively  (Kim  et  al., 
1993). 

The  PgILg  spectral  amplitude  ratio  using  vertical-component  records  in  a  broad  high 
frequency  band  of  4-18  Hz  provides  discrimination  power  with  a  total  misclassification  probability 
of  about  6.3-7. 1%.  In  a  lower  frequency  band,  2-10  Hz,  PgILg  ratios  show  poor  results  and  the 
separation  between  explosions  and  earthquakes  is  less  clear  with  a  total  misclassification 
probability  of  about  14.3%  (Table  4).  However,  when  a  single  PgILg  ratio  at  12  Hz  is  added  in 
the  low  frequency  band,  the  new  frequency  band,  2-12  Hz,  improves  the  discrimination  power  by 
up  to  about  6%  for  vertical-  as  well  as  three-component  data  (Table  4).  This  significant 
improvement,  which  is  due  to  the  addition  of  a  single  12  Hz  PgILg  ratio,  illustrates  the 
importance  of  high-frequency  data  in  regional  seismic  discrimination. 

The  PgILg  spectral  ratios  using  three-component  regional  records,  corrected  for  the  free 
surface  effect,  improve  the  discrimination  power  of  the  spectral  ratio  method.  In  the  frequency 
band  8—18  Hz,  the  free  surface  corrected,  three-component  PgILg  spectral  ratio  provides 
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discrimination  power  with  a  total  misclassification  probability  of  only  2.6%. 

In  the  case  that  free  surface  effects  cannot  be  corrected  with  confidence,  then  the  PgILg 
spectral  ratios  of  rotated,  three-component,  regional  records  may  be  used  as  a  discriminant.  The 
use  of  rotated,  three-component  records  improves  the  discrimination  power  over  the  vertical- 
component.  In  the  frequency  band  of  8-18  Hz,  the  three-component  PgILg  spectral  ratio 
provides  discrimination  power  with  a  total  misclassification  probability  of  about  3%. 

The  area  around  the  Caucasus  Network  was  not  an  ideal  region  to  make  free  surface 
corrections,  because  of  site  effects  associated  with  low  velocity  near-surface  layering.  Application 
of  the  free  surface  correction  to  data  acquired  at  a  good  hard  rock  site  gives  a  better  separation  of 
P,  SV  and  SH  signals  on  different  traces.  For  such  sites,  we  would  therefore  expect  an  even 
better  discrimination  performance  for  the  free  surface  corrected  Pg/Lg  ratio. 

Preliminary  analysis  of  applying  distance  corrections  to  the  vertical-component  Pg/Lg 
spectral  ratio  indicates  that  use  of  the  distance  corrected  spectral  ratios  improves  the  discrimination 
power  of  high-frequency  regional  data  by  about  4%  over  the  use  of  vertical-component  ratios 
(92.9%  versus  97.0%).  However,  accurate  information  regarding  the  geometrical  spreading  and 
anelastic  attenuation  of  the  regional  phases,  such  as  Pg  and  Lg  at  high  frequency,  are  required  for 
applying  distance  corrections  with  confidence  for  the  region  in  question.  Lacking  more  detailed 
information  on  amplitude  decay  of  regional  phases  with  distance,  in  particular  frequency  dependent 
Q  factors  for  P  (or  Pg)  and  Lg  waves,  in  the  northern  Caucasus  region,  we  applied  distance 
corrections  to  only  the  vertical-component  Pg/Lg  ratios  in  this  study. 

We  have  shown  that  the  Pg/Lg  ratio  is  an  adequate  discriminant  for  explosions  from 
earthquakes  (magnitude  smaller  than  4.5)  in  southern  Russia,  and  that  the  discrimination  is 
improved  by  the  use  of  three-component  data  together  with  a  free  surface  correction.  The  stability 
of  the  Pg/Lg  ratio  in  discriminating  multiple-hole,  ripple-fired  explosions  from  small  regional 
earthquakes  is  significant  for  seismic  monitoring  in  southern  Russia.  Our  results  demonstrate  the 
importance  of  data  at  frequencies  up  to  at  least  20  Hz. 

One  basis  for  the  empirical  success  of  our  proposed  high-frequency  Pg/Lg  ratio  method  is 
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that,  if  selected  spectral  amplitude  reinforcement  occurs  at  the  source  (as  is  the  case  for  ripple-fired 
blasts),  it  will  affect  both  early  P  phases  as  well  as  later-arriving  Lg  phases.  Such  spectral 
scalloping  will  be  largely  canceled  in  the  ratio  we  have  used.  The  ratio  is  a  robust  discriminant  to 
the  extent  that  there  are  cancellation  of  other  characteristics,  such  as  event  size,  comer  frequency, 
some  effects  of  focal  depth,  and  instrument  response. 

The  problem  of  certain  earthquake  radiation  patterns  is  also  minimized  by  using  the  network 
averaged  ratios  and  three-component  records,  for  example,  a  strong  Pg  excitation  but  weak  Lg 
excitation  from  a  near  vertical  strike-slip  earthquake  (see  e.g.,  Bouchon,  1981;  Kim,  1987)  which 
may  diminish  the  discrimination  power  of  the  PgILg  ratio  method  at  some  station  azimuths,  as 
indicated  by  Lilwall  (1988). 

A  key  to  the  portability  of  the  PgILg  spectral  ratio  discriminant  to  various  parts  of  the  world 
is  the  correction  of  regional  signals  for  their  source-receiver  path  effects.  We  have  shown  that  the 
free  surface  effect  can  be  easily  removed  from  the  regional  records  and  that  free  surface  corrected, 
PgILg  ratios  using  three-component  records  improve  the  discrimination  power  of  the  PgILg 
ratio  method.  Further,  the  dependence  of  the  frequency  content  of  P  and  S  waves  on  specific 
propagation  paths  and  local  structure  at  the  receiver  is  reduced  by  using  the  network  averaged 
PgILg  ratio  (see  also  Kim  et  al.,  1993). 

The  PgILg  spectral  ratios  may  be  taken  for  different  frequencies,  for  instance  Pg(high 
frequency )/Lg(low  frequency)  ratios  or  ratios  at  different  frequencies  within  the  same  phase  as 
Lg(high  frequency )/Lg(low  frequency)  ratio.  Such  cross  spectral  measures  may  be  attractive  for 
discriminating  between  distant  earthquakes  and  explosions,  from  a  few  hundred  to  a  couple  of 
thousand  kilometers.  However  we  need  to  correct  for  instrument  responses;  and  propagation  path 
effects  -  notably  the  frequency  dependent  anelastic  attenuation  -  must  be  taken  into  account  for  the 
cross  spectral  ratios  to  be  meaningful.  The  multivariate  analysis  employed  in  this  study  implicitly 
takes  such  a  cross  spectral  measure  when  a  broad  enough  frequency  band  is  chosen  for  the  input 
variables. 

It  has  been  reported  by  some  researchers  (e.g.,  Fisk,  1993;  Taylor,  1996)  as  well  as  in  this 
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study  that  the  distance  corrections  -  primarily  for  anelastic  attenuation  along  the  paths  -  applied  for 
both  the  P  and  S  signals  enhance  the  discriminant  power  of  PgILg  ratios.  However,  we  believe 
that  it  is  difficult  to  correct  signals  for  the  effects  of  distance  with  confidence,  unless  special  studies 
are  carried  out,  for  example  with  empirical  Green’s  functions.  The  PgILg  spectral  amplitude  ratio 
method  we  propose  is  sufficiently  reliable  and  robust  that  it  can  be  used  in  discriminating  chemical 
explosions  (especially  numerous  mining  and  quarry  blasts)  from  small  regional  earthquakes  in  the 
routine  analysis  of  data  from  regional  earthquake  monitoring  networks. 

In  the  case  of  regions  with  high  attenuation  such  as  Nevada,  the  PgILg  spectral  ratio 
method  may  be  applicable  with  somewhat  diminished  discrimination  power  when  compared  with 
stable  regions,  since  the  frequency  band  available  for  regional  signals  in  such  a  region  may  be 
limited  to  up  to  about  10  Hz.  A  comparison  of  performance  between  low-  and  high-frequency 
band  PgILg  spectral  ratios  summarized  in  Table  4  indicates  that  the  method  should  work  in  such 
regions  (using  frequencies  up  to  10  Hz)  with  a  total  misclassification  probability  of  about  10-15%. 
It  may  be  emphasized  that  although  we  have  presented  results  of  frequency  bands  with  the  best 
discrimination  power  for  four  different  cases,  any  combination  of  frequency  bands  across  the 
range  2  -  18  Hz  provides  fairly  good  results  (Table  4).  For  instance,  for  vertical-component  data, 
any  six  discrete  frequencies  between  2-18  Hz  provide  a  total  misclassification  probability  of 
between  7.1  and  8.4%  (Table  4). 
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Figure  Captions 

Fig.  1.  Locations  of  the  Caucasus  Network  stations  (solid  squares)  and  near  regional  (A  <  300 
km)  seismic  events  {shaded  circles)  located  by  the  Caucasus  Network  during  1992  are  plotted  on 
the  topographic  relief  map.  Notice  the  clusters  of  epicenters  near  43.4°N:42.8°E  (Tyrnauz  mine; 
Table  1),  44°N:42°E  (Ust-Djeguta  quarry;  Table  1)  and  42.5°N;43.5°E  (epicentral  area  of  04/29/91 
Ms  =  1  Georgian  earthquake).  Symbol  sizes  are  proportional  to  the  magnitudes  of  the  events. 
Three  major  tectonic  features  in  the  region  -  Russian  Platform,  Greater  Caucasus  and  Lesser 
Caucasus  -  are  indicated. 

Fig.  2.  A  comparison  of  frequency  contents  of  the  regional  seismic  signals  from  an  earthquake 
and  a  ripple-fired  quarry  blast,  (a)  spectrogram  of  the  earthquake  record  shows  strong  S  wave 
amplitude  at  low  frequencies  (below  about  5  Hz),  and  at  higher  frequencies  (10-15  Hz),  P  and  S 
waves  have  comparable  amplitude,  (b)  spectrogram  of  the  quarry  blast  record  shows  predominant 
P  wave  energy  over  a  wide  frequency  band  3-20  Hz,  but  weak  S  waves  in  this  band.  The 
spectrogram  of  earthquake  record  shows  random  distribution  of  seismic  wave  energy  while  the 
explosion  signal  shows  clear  spectral  banding  due  to  source  multiplicity  (ripple-firing;  a  spectral 
band  at  about  22  Hz  is  due  to  noise). 

Fig.  3.  Locations  of  earthquakes  (circles),  mines  and  quarries  identified  and  used  in  this  study 
(solid  triangles),  other  known  quarries  (open  triangles)  and  Caucasus  Network  stations  (solid 
squares).  Shaded  region  represents  areas  over  2,000  m  elevation  topped  by  El’brus  Mt.  (5633  m) 
which  is  the  highest  peak  in  Europe  west  of  the  Urals.  Circle  size  is  proportional  to  the  magnitude 
of  the  earthquakes.  Numbers  for  the  quarries  and  earthquakes  refer  to  the  event  ids  given  in  Table 
1  and  Table  2. 


Fig.  4.  Network  averaged  vertical-component  log^Q(Pg/Lg)  spectral  amplitude  ratios  at  12 
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discrete  frequency  points  used  in  discrimination  analysis  are  plotted  for  earthquakes  {shaded 
circles)  and  explosions  {open  triangles).  A  mean  value  of  each  discrete  frequency  point  is  plotted 
for  earthquakes  {open  circles)  and  explosions  {solid  triangles)  with  their  arms  representing  the 
scatter  (s.d.). 

Fig.  5.  Discriminant  scores  of  earthquakes  {shaded  circles)  and  explosions  {open  triangles)  of 
the  sample  data  are  plotted  with  their  mean  network  averaged  logjQ(Pg/Lg)  ratios.  Note  that  the 
two  populations  are  also  separated  by  mean  vertical-component  logjQ(Pg/Lg)  ratio  «  0.25. 
Vertical  lines  denoted  as  Ex  and  Eq  are  the  projection  of  the  multivariate  means  of  the  earthquake 
and  explosion  population,  respectively.  The  vertical  line  Dq  is  the  classification  line.  Two 
earthquakes  and  an  explosion  are  incorrectly  classified  and  the  total  misclassification  probability  is 
7.1  %. 

Fig.  6.  Network  averaged  three-component  logjQ(Pg/Lg)  spectral  amplitude  ratios  at  12  discrete 
frequency  points  used  in  discrimination  analysis  are  plotted  for  earthquakes  {shaded  circles)  and 
explosions  {open  triangles).  Other  notation  is  the  same  as  in  Figure  4. 

Fig.  7.  Discriminant  scores  of  earthquakes  {shaded  circles)  and  explosions  {open  triangles)  of 
the  sample  data  are  plotted  with  their  mean  network  averaged  logjQ(Pg/Lg)  ratios.  All 
earthquakes  and  explosions  are  correctly  classified  and  the  total  misclassification  probability  is  3.0 
%.  Each  event  is  indicated  by  its  event  id  numbers  given  in  Table  2  &  3  and  in  Figure  3.  Notice 
that  distant  earthquakes,  #2,  #8  and  #16,  have  relatively  high  spectral  ratios.  Other  notation  is  the 
same  as  in  Figure  5. 

Fig.  8.  Sketch  of  ground  motions  in  a  homogeneous  half-space  for  a  P-wave  incident  upon  a  free 
surface  with  unit  amplitude  and  an  angle  of  incidence  i  «  33.37°  (=  0.1  s/km).  Surface  P  and  S 
velocities  are:  oc  =  5.5  km/s  and  =  3.1  km/s,  respectively.  The  total  ground  motion  of  the  free 
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surface  {shaded  line)  itself  is  shown  above  the  surface  level  z  =  0  and  has  amplitude  2.03  times 
the  incident  wave  with  an  apparent  angle  of  incidence  i'  =  36.12°.  Scattered  waves  PP  and  PS 
which  are  downgoing  waves,  are  plotted  above  the  level  z  =  0  indicating  the  directions  of  true 
particle  motion.  The  amplitudes  of  particle  motions  are  indicated  by  their  length. 

Fig.  9.  {Upper  three  traces)  Rotated  three-component  records  from  an  explosion  in  Tymauz  mine 
(Z,  R  and  T).  The  P  waves  are  on  the  Z-component,  while  S  waves  are  dominantly  on  the  R 
and  T  components.  {Lower  three  traces)  Composite  incident  wave  vector  component  traces  {P, 
SV  and  SH)  produced  by  applying  free  surface  correction  operators  successively  along  the  trace. 
The  breaks  between  different  slowness  bands  are  marked  by  vertical  bars  and  the  slowness  (in 
s/km)  used  in  each  group  velocity  window  is  indicated.  All  traces  are  plotted  with  the  same 
amplitude  scale.  Trace  peak  to  peak  amplitudes  in  p/s  are  given  in  each  trace. 

Fig.  10.  Network  averaged  log^Q(Fg/Lg)  spectral  amplitude  ratios  at  12  discrete  frequency 
points,  from  the  free  surface  corrected  three-component  records  (P-,  SV-,  and  5//-seismogram) 
used  in  discrimination  analysis,  are  plotted  for  earthquakes  {shaded  circles)  and  explosions  {open 
triangles).  Other  notation  is  the  same  as  in  Figure  4. 

Fig.  11.  Discriminant  scores  of  earthquakes  {shaded  circles)  and  explosions  {triangles),  of  the 
sample  data  from  the  free  surface  corrected  three-component  records,  are  plotted  with  their  mean 
network  averaged  log|Q(Pg/Lg)  ratios.  Other  notation  is  the  same  as  in  Figure  5. 

Fig.  12.  The  vertical-component  mean  high-frequency  (8-18  Hz)  PgILg  ratios  of  82  earthquake 
records  {open  circles)  and  78  explosion  Tecords{open  triangles)  are  plotted  against  distance. 
Corresponding  distance  corrected  spectral  ratios  are  plotted  with  shaded  symbols.  An  average 
trend,  calculated  from  each  of  the  20  km  distance  bins,  is  superposed  {thick  lines)  and  an  average 
trend  after  the  distance  correction  is  plotted  for  comparison  {thin  lines).  The  mean  ratio  in  each 
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distance  bin  is  plotted  at  the  mean  distance  in  each  bin  for  earthquake  {open  circle)  and  for 
explosion  (open  triangle)  and  their  sizes  are  proportional  to  the  density  of  the  bin. 

Fig.  13.  Distance  correction  factors  for  logjQ(Fg/Lg)  spectral  ratios  for  the  northern  Caucasus 
region  are  plotted  as  a  function  of  frequencies  and  distances.  Correction  factors  for  various  paths 
are  indicated:  shield  path  {dotted  lines),  mountain  {dashed  lines)  and  mixed  path  -  path  crossing 
both  regions  {solid  lines).  The  correction  factor  includes  a  differential  geometrical  spreading  of 
AO.5  (for  A  >  100  km)  and  attenuation  coefficient  (Yp  -  y^). 

Fig.  14.  Network  averaged  vertical-component  logjo(Pg/Lg)  spectral  ratios  corrected  for  the 
distance  are  plotted  for  earthquakes  {shaded  circles)  and  explosions  {open  triangles).  At  each 
discrete  frequency,  mean  values  of  distance  uncorrected  ratios  are  also  plotted  for  earthquakes 
{open  circles)  mdi  explosions  {solid  triangles)  for  comparison.  Other  notation  is  the  same  as  in 
Figure  4. 

Fig.  15.  Discriminant  scores  of  earthquakes  {shaded  circles)  and  explosions  {open  triangles)  of 
the  distance  corrected  vertical-component  data  are  plotted  with  their  mean  network  averaged 
^ogio(Fg/Lg)  ratios.  Note  that  these  two  populations  are  also  well  separated  by  mean 
^oSio{P8lLg)  ratios.  One  earthquake  is  incorrectly  classified  and  the  total  misclassification 
probability  is  3.0%.  Each  event  is  identified  by  its  id  number  given  in  Table  2  &  3.  Other  notation 
is  the  same  as  in  Figure  5. 
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Table  1.  Quarries  and  Mines  around  Kislovodsk,  Russia^^^ 


Id 

Quarry 

Lat. 

Long. 

Id 

Quarry 

Lat. 

Long. 

# 

name 

(°N) 

(°E) 

# 

name 

(°N) 

(°E) 

1 

Alenovsky 

44.20 

42.31 

8 

T  sementny-Zavod 

44.11 

42.04 

2 

Andezit 

43.74 

41.91 

9 

Tymauzt 

43.39 

42.87 

3 

Byely-Mramoz 

43.48 

41.76 

10 

Ust-Djeguta 

44.09 

42.02 

4 

Cher-Mramoz 

43.56 

41.79 

11 

Verblud 

44.18 

42.96 

5 

Zhakot 

43.94 

41.83 

12 

Zamok 

43.89 

42.62 

6 

Grznimy 

43.35 

42.19 

13 

Zmeika 

44.18 

43.11 

7 

Kinjal 

44.27 

43.01 

14 

Unidentified 

43.50 

43.91 

Based  on  local  draft  map  (Dr.  Vitaly  Khalturin,  person,  comm.,  1994); 
t  F.  Riviere-Barbier,  Center  for  Seismic  Studies,  Washington  D.  C.  (person,  comm.,  March 
1994). 
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Table  2.  List  of  earthquakes  used  in  this  study 


Id 

# 

Origin  time 

Lat. 

Long. 

Depth 

mb 

Distance 

3-component 

yrmody  hh:mm:sec 

CN) 

(°E) 

(km) 

JSPCt 

(km) 

Score 

Pg/Lg 

N 

1 

92  02  16 

13:43:14.8 

44.398 

42.438 

16 

0.5 

40.7 

11.10 

-0.11 

2 

2 

92  02  20 

20:18:55.1 

42.526 

43.280 

3.4 

197.2 

5.34 

0.12 

3 

3 

92  02  26 

06:13:28.9 

44.317 

42.536 

18 

-0.2 

53.9 

13.38 

-0.43 

2 

4 

92  02  28 

01:16:53.7 

44.394 

42.632 

18 

0.9 

50.3 

11.50 

-0.18 

3 

5 

92  03  07 

03:06:42.9 

44.069 

43.024 

30 

- 

43.6 

10.31 

-0.23 

4 

6 

92  03  07 

03:08:14.0 

44.075 

43.028 

11 

- 

43.6 

9.95 

-0.21 

4 

7 

92  03  22 

05:10:00.9 

43.279 

41.411 

2.9 

144.6 

2.27 

0.05 

3 

8 

92  03  27 

19:21:04.8 

42.466 

43.725 

33 

4.8tt 

217.7 

0.32 

0.27 

3 

9 

92  03  28 

09:39:04.8 

44.089 

42.770 

12 

- 

48.8 

9.42 

-0.24 

4 

10 

92  04  07 

20:31:09.7 

43.206 

43.491 

15 

4.5 

134.5 

5.18 

-0.17 

2 

11 

92  04  09 

17:34:00.9 

42.594 

43.361 

2.8 

175.0 

2.88 

-0.01 

3 

12 

92  04  12 

02:16:30.9 

44.250 

42.617 

15 

0.0 

45.2 

6.70 

-0.20 

3 

13 

92  04  12 

20:00:56.2 

43.923 

43.081 

28 

-0.1 

43.1 

13.11 

-0.36 

3 

14 

92  04  16 

14:18:57.5 

43.727 

43.035 

6 

- 

60.4 

10.80 

-0.45 

3 

15 

92  05  09 

20:45:06.8 

42.670 

43.519 

2.9 

173.6 

8.50 

-0.09 

3 

16 

92  06  21 

21:14:24.6 

43.577 

44.625 

31 

2.8 

129.3 

2.38 

0.16 

2 

17 

92  07  05 

17:34:07.0 

45.298 

42.212 

13 

3.8 

149.8 

6.01 

-0.03 

3 

18 

92  07  11 

14:39:17.8 

42.563 

43.505 

2.6 

184.8 

8.47 

-0.10 

3 

19 

92  07  16 

17:15:33.3 

43.433 

43.194 

2.9 

91.2 

0.47 

0.01 

3 

20 

92  07  16 

17:20:09.6 

43.440 

43.141 

98.9 

2.43 

-0.02 

2 

21 

92  07  17 

19:42:03.7 

43.423 

43.166 

8 

77.5 

2.95 

-0.09 

2 

33 


22 

92  07  23 

01:07:02.8 

44.164 

42.609 

10 

37.1 

3.59 

-0.19 

3 

23 

92  07  23 

01:13:15.0 

44.104 

42.890 

17 

0.1 

42.5 

13.04 

-0.51 

2 

24 

92  08  12 

13:23:18.6 

43.440 

43.578 

22 

3.3 

94.6 

6.75 

-0.13 

4 

25 

92  08  21 

20:12:46.7 

44.326 

41.757 

1.9 

99.4 

9.07 

-0.32 

3 

^  Magnitude  given  in  CAUCASUS  NETWORK  Information  Product,  Technical  Reference 
Manual,  IRIS/Joint  Seismic  Program  Center  (JSPC),  Univ.  Colorado  and  Lamont-Doherty  Earth 
Observatory,  April  1994;  tt  Magnitude  given  in  PDE  monthly  list;  Distance=mean  epicentral 
distance  for  each  event  used  in  3-component  analysis;  Score=discriminant  score  from  3-component 
analysis  (see  Figure  7);  Pg/Lg=mean  spectral  ratio  used  in  3-component  analysis;  N=number  of 
records  used  in  3-component  analysis. 


Table  3.  List  of  quarry  blasts  used  in  this  study 


Id  Origin  time  /nj.  Distance  3-component  analysis 


# 

yr  mo  dy 

hh:mm:sec 

JSPC 

(km) 

Score 

Pg/Lg 

N 

Tymauz  mine  (43.386°N,  42.87 1°E) 

1 

92  02  23 

07:02:43.4 

3.3 

89.5 

-8.07 

0.53 

4 

2 

92  03  11 

13:39:01.5 

3.2 

96.8 

-15.10 

0.64 

3 

3 

92  03  22 

07:29  38.7 

86.2 

-2.09 

0.34 

2 

4 

92  03  29 

06:42:46.5 

2.7 

86.2 

-6.56 

0.40 

2 

5 

92  04  02 

12:19:39.3 

80.0 

-10.79 

0.51 

3 

6 

92  04  09 

13:02:38.9 

80.0 

-4.79 

037 

3 

7 

92  04  12t 

08:04:34.1 

80.0 

-7.66 

0.44 

3 

8 

92  04  19 

03:53:37.0 

79.7 

-8.53 

0.46 

3 

9 

92  04  23 

13:28:42.4 

2.0 

86.2 

-5.12 

0.44 

2 

10 

92  06  07 

05:39:56.2 

1.5 

65.6 

-7.02 

0.35 

2 

11 

92  06  07 

05:55:40.7 

2.1 

83.0 

-10.14 

0.38 

3 

12 

92  06  14 

04:37:18.3 

83.0 

-7.81 

0.44 

3 

13 

92  07  05 

06:35:03.4 

65.6 

-5.78 

0.36 

2 

14 

92  07  12 

06:35:54.8 

66.0 

-5.50 

0.20 

2 

15(a)  93  08  22t 

05:09:38.4 

18.9tt 

53.4 

-11.41 

0.35 

3 

16(b)  93  08  22t 

06:24:24.4 

24.0tt 

53.4 

-6.70 

0.27 

3 

17(c)  93  08  29t 

05:23:58.6 

10.2tt 

53.4 

-12.52 

0.29 

3 

Ust  Djeguta  (44. 10°N,  42.05°E) 

56.2  -7.66  0.23  1 


18  92  01  27  09:32:45.1 
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19 

92  02  19 

11:34:41.7 

52.2 

-6.74 

0.26 

3 

20 

92  03  16 

11:37:17.1 

71.9 

-6.15 

0.24 

1 

21 

92  04  14 

11:02:09.7 

69.8 

-1.67 

0.13 

3 

22 

92  04  21 

11:15:58.2 

69.8 

-2.32 

0.08 

3 

23 

92  06  09 

14:03:37.8 

69.3 

-8.40 

0.33 

1 

24 

92  06  24 

09:41:05.9 

85.3 

-6.43 

0.22 

2 

Zhako-Krasnogorskaya  (43.941°N,  41 

.831°E) 

25 

92  05  08t 

11:23:50.0 

63.7 

-0.99 

0.08 

2 

^  Locations  and  origin  times  provided  by  F.  Riviere-Barbier,  Center  for  Seismic 
Studies,  Washington  D.  C.  (person,  comm.,  March  1994);  tt  metric  tons  (B. 
Stump,  person,  comm..  May  1994);  (a)  underground  1-3  shots;  (b)  surface 
ripple-fired;  (c)  underground  1-3  shots;  Distance=mean  distance  for  each  event 
used  in  3-component  analysis.  Other  notation  is  the  same  as  in  Table  2. 
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Table  4.  Comparison  between  different  frequency  bands  in  application  of  the  PgILg  discriminant. 


Case 

2  - 10  Hz 

10- 18  Hz 

2- 

12  Hz 

8  - 18  Hz 

Pt 

A2tt 

P 

A2 

P 

A2 

P 

A2 

(%) 

(%) 

(%) 

(%) 

Vertical-component 

85.70 

4.5 

92.10 

8.0 

91.56 

7.6 

92.95 

8.7 

Three-component 

87.96 

5.5 

96.96 

14.1 

90.29 

6.7 

96.96 

14.1 

Free  surface  corrected 

89.85 

6.5 

97.14 

14.5 

93.01 

8.7 

97.42 

15.1 

Distance  corrected 

89.02 

6.0 

96.01 

12.3 

96.50 

13.1 

96.99 

14.1 

t  probability;  tt  Mahalanobis  D-squared  measure. 


Network  Averaged  LoQ^^{Pg/Lg),  Z-component 


Network  Mean  Log  (P^/Lg),  Z-componen 


Network  Averaged  Log^Q(P^/Lg),  3-component 


03/29/92  06:42:46.5,  43.386"N,  42.871  "E,  m^=  2.7,  Tyrnauz  mine 


.08 1  .14 


.25 


.34 


MIC 

P 

2.18  p,/s 


.37  s/km 


108.0  km 
az=  337.3° 
baz=156.9° 


SH 

1.75 


"#lhVl(||WY^Y'AAA/VvAM/V^^ 


[iiii|iiii|iiii|  mi  1 1  n  iTTn-'i' I  I  i  i  i  |  i  i  i  i — | — i — i — i — i - 1 - 1 - 1 - 1 - 1 - 1 

7  6  5  4  3  2 

Group  Velocity  (km/sec) 


Network  Averaged  LoQ^^{Pg/Lg),  Free  surface  corr 


Mean  Log^Q(Pg/L^),  8-18  Hz,  Z-component 


A 


0.84 


0.6 

0.4 

0.2 


0.0  H 


-0.2 


-0.4- 

-0.6- 


-0.8 


0 


b 
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Network  Averaged  LoQ.APg/Lg),  Distance  corrected 


